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Abstract 

We numerically investigate density perturbations generated in the smooth hybrid new inflation 
model, a kind of double inflation model that is designed to reproduce the running spectral index 
suggested by the WMAP results. We confirm that this model provides the running spectral index 
within la range of the three year WMAP result. In addition, we find a sharp and strong peak on the 
spectrum of primordial curvature perturbation at small scales. This originates from amplification 
of fluctuation in the first inflaton fields due to parametric resonance, which takes place in the 
oscillatory phase between two inflationary regime. Formation probability of primordial black holes 
(PBHs) is discussed as a consequence of such peak. 
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I. INTRODUCTION 



The observation of Wilkinson Microwave Anisotropy Probe (WMAP) has successfully 
determined the cosmological parameters and nature of density fluctuations with high preci- 
sion. The WMAP result is quite consistent with the prediction of the inflationary cosmology, 
i.e., the flat universe with almost scale-invariant adiabatic density fluctuation. If one take 
a closer look at the shape of the power spectrum of the density fluctuations obtained by 
WMAP, however, small deviation from scale invariance is found. In fact, the WMAP 1st 
year result suggested that the data fits to the power spectrum with running spectral index 
n s , although the statistical significance was not high enough [3]. There has been a 
renewed interest in the running spectral index because it is favored by the recently published 



three year WMAP result 



4j which gives dn s /d\nk = — 0.055^'^- In Ref. no evidence 



of running was found from the combining analysis of the 1st year WMAP and the improved 
data of Ly-a forests. However, as mentioned in 4|, the three year WMAP result is not in 
good agreement with this Ly-a data. Then, it seems premature to adopt the power spec- 
trum obtained from Ly-a forests. Therefore it is worth studying the inflation model which 
provides the running spectral index. 

It is not an easy task to build an inflation model which produces density perturbations 



A 



ter the release of the 1st year WMAP 



9J, smooth hybrid+new [10|]) were pro- 



whose power spectrum has a running index {(], 7 
data, the double inflation models (hybrid+new 
posed. However, since the power spectrum was obtained by analytical method in those 
works, the precise form of the spectrum on scales corresponding to transition from one infla- 
tion to another was not clear, which makes it difficult to compare the theoretical prediction 
with observations. 

In this paper, therefore, we consider a double inflation model and calculate the power 
spectrum by numerical integration of evolution equations for density fluctuations. We adopt 
the smooth hybrid +new inflation clS Sb double inflation model. The model consists 
of smooth hybrid inflation [ll| and new inflation [l2, Q], both of which are based on 
supergravity. The running index is realized by the smooth hybrid inflation. The new inflation 
is necessary because the hybrid inflation with a large running index only has small e-folds 
iV ~ 10 and the remaining e-folds for successful inflation are provided by the new inflation. 
We show that the model can give an appropriate power spectrum which are consistent with 
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WMAP result. Moreover, we find that large density fluctuations are produced through 



parametric resonance 



15 



Il6l . Il7j of the inflaton fields at the transition epoch, which leads 
to a sharp peak on small scales in the spectrum and formation of primordial black holes 
(PBHs). 

This paper is organized as follows: In Sec. II, we introduce the smooth hybrid new in- 
flation. The results of our numerical calculation on the power spectrum of the density 
perturbation are shown in Sec. III. In Sec. IV, we investigate the parametric resonance in this 
inflationary model which results in a sharp peak found in Sec. III. In Sec.V, PBH formation 
is discussed as a consequence of a resonant peak. Sec. VI is devoted to a discussion. 

In this paper, we set the reduced Planck scale M G = 2.4 x 10 18 GeV to be unity unless 
otherwise stated. 



II. SMOOTH HYBRID NEW INFLATION MODEL 



In this section, we briefly review the smooth hybrid new inflation model lOj, which is 
based on supergravity. Its superpotential is given by 

W = W H + W N , (1) 

where Wh (Wn) is the superpotential responsible for smooth hybrid (new) inflation. Wh is 
written as 

( (W) m A 

w * = s [-f + h^)- (2) 

Here, $ and ^ are a conjugate pair of superfields transforming as nontrivial representations 
of some gauge group G. S is a superfield whose scalar component is the inflaton and 
transforms as a singlet under G. Moreover, Wh has two symmetries: one is an i?-symmetry 
under which S — > e 2m ,\l/ — > \l/ and ^ — ► and the other is a discrete Z m symmetry 
under which the combination tyty has unit charge. M sets a cutoff scale which controls 
nonrenormalizable terms in Wh- We include possible coupling constants in the definition of 
M. n sets the scale of the smooth hybrid inflation. 
The superpotential Wn is given by 

W N = v 2 ® — $ n+1 . (3) 

n+1 



<3> is the inflaton superfield for the new inflation and has a discrete i?-symmetry, Z2 U r [13 ]. 
v , which satisfies v <C fi, is the scale of the new inflation, and g is a coupling constant of 
nonrenormalizable terms in W^. 
The Kahler potential is given by 

K = K H + K N , (4) 



K N 



\S\ 2 + |^| 2 + |^| 
11 + 4 1 1 ' 



(5) 
(6) 



where Cat is a constant smaller than unity. 

From W and K , we can derive the scalar potential. We assume .D-flatness, which coincides 
with the steepest descent direction in the F-term contribution, and consider only F-term 
contribution. The scalar potential is given by 



V = exp 



151 



|tf|2 + |vj>| 2 + |$| 2 + ^|$| 4 



X 



:i+\s\ 2 ) v+ 



M 2 ( m_1 ) 



n + 1 



+ 
+ 
+ 



m—l 



n + 1 



5'(^^) m _ , 

mW M 2(m-1) ' ^ 



2 , W 



M2{m~l) 



1 



l + Cjvl^] 2 



l$l 
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n + 1 
Cat|$| 4 



V 2( 1+ | $ |2 + ^V| $ |4 j _ (1 

2 1 1 / V n+1 2 n+1 



( 1 + ^|$| 2 



-/I 



M 2 ( m_1 ) 



3IWI 



(7) 



Here, the scalar components of the superfields are denoted by the same symbols as the 
corresponding superfields. Performing adequate transformations allowed by the symmetries, 
the complex scalar fields are changed into real scalar fields: 



o = V2ReS, $ = 2Re^ = 2Re^, = v^Re^. 



(8) 



In what follows, we will use these real scalar fields. 
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A. Smooth hybrid inflation 



First, we make the assumption that initially |<r| is sufficiently large though \a\ < 1 is 
satisfied, and that ip and <fi are set around local minimum. Indeed, if \a\ is sufficiently large, 
ip and have effective masses larger than the Hubble parameter H, so that they roll down to 
their respective minima quickly. Thus, the effective potential of a determines not only the 
dynamics of the smooth hybrid inflation but also primordial density fluctuations generated 



during the smooth hybrid inflation 18|. Retaining only relevant terms, it is given by 



, , 2 M 2 u 2 a 4 . 
V HeS (a) ~ ( 1 - -— f - + _ + -..). (9) 



Hereafter, we consider only the case with m = 2. As long as y/Mjx <C |er| -C 1, the effective 
potential is dominated by the false vacuum energy /i 4 , and hence inflation takes place. The 
Hubble parameter H is almost constant H ~ Hh = /i 2 /"\/3- The second term of Eq.([9]), 
which originates from nonrenormalizable term of superpotential, has a negative curvature. 
On the other hand, the third term of Eq.Q, which comes from supergravity correction, has 
a positive curvature. Then, for modes crossing the horizon while the third term dominates 
the dynamics, the spectrum of curvature perturbation has a spectral index n s > 1. In 
opposition, for modes crossing the horizon while the second term dominates, the spectrum 
of curvature perturbation has a spectral index n s < 1. In this way, this smooth hybrid 
inflation model can generate the spectrum with the running spectral index suggested by 
WMAP results. 

We estimate the amplitude of curvature perturbation TZ, spectral index n s and its running 



dn s /d\nk, up to second order of slow-roll parameters (Ref. 1Q(|). According to WMAP three 



year result, these cosmo logical parameters are constrained by 

U 2 (k ) = 24+i x 10 -h> ns{ko) = i.050+^ 4 , ^ = -0.055+™! (10) 
at la level, where k = 0.002 [Mpc -1 ]. We can determine the parameters requiring that 



the spectrum satisfies these constraints. However, as shown in Ref. [10], the number of 
e-foldings Nh after the mode with k = 0.002 [Mpc -1 ] crosses the horizon is estimated to be 
N[{ < 10. This is too small to solve the horizon and the flatness problems, and hence the 
second inflation is needed. 

In addition, we must require that Nh is sufficiently large: Nh > 10, for, as we will see 
later, the perturbations produced by the second inflation has a very large amplitude. The 



scale of those perturbations should be sufficiently small, say k > l[Mpc _1 ], in order not to 
conflict with observations. It is difficult to meet this request while reproducing the best-fit 
values of WMAP three year result, because the latter requires large running. If we allow la 
range of WMAP three year constraint, we can reproduce such a spectrum. We choose the 
following parameters: 

H = 2.04 x 10" 3 , M = 1.17. (11) 
This reproduces best-fit values of TZ and n s , but dn s /d\nk within the la range. 

H(ko) = 4.9 x 10" 5 , n s {k ) = 1.053, —\ = -0.032 (12) 

& 111 K 

at k — 0.002 [Mpc -1 ], as the result of numerical calculation described in Sec. III. 

Note that the present model produces negligibly small tensor modes. In fact, the tensor 
to scalar ratio r is estimated by a slow-roll parameter as r = lQe H < 10~ 3 . 



B. Oscillatory phase 

After the smooth hybrid inflation, the oscillatory regime sets in. a and ip oscillate around 
their respective minima, cr min = and ^ min = l^fJiM . Hereafter, we replace tp — * ip m m + i ) - 
Around these minima, their effective masses m a and are given by 



^ = Vi7' ^ = v^ +16/i4 - (13) 

These scalar fields a and tp undergo damped oscillations, whose amplitude decreases asymp- 
totically in proportion to t^ 1 oc a~ 3//2 . The total energy of oscillating fields a and ip decreases 
as a~ 3 , like non-relativistic matter. Eventually, contribution of 0, which is about t> 4 , dom- 
inates the total energy, and the new inflation starts. The duration of oscillatory phase in 
terms of the scale factor can be estimated as 

l n ^i~!l n ^. (14) 
a c 3 v 

The subscript "c" indicates that the value is evaluated at the end of smooth hybrid inflation, 
and the subscript "ini" indicates that the value is evaluated at the beginning of the new 
inflation. 
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On the other hand, due to interaction terms between a, ip and 0, a ls° has an effective 
mass rriff). Averaged for a time-scale sufficiently longer than the period of oscillation of a 
and ip, one can estimate m^ as 

ml ~ -if 2 . (15) 

The effective mass of is larger than H, so it shows oscillatory behavior with a very long 
period Q. Toe a mplitU de of th- osci llat io n deca y8 in proportion to tx a 3 ' 4 . At the 
beginning of this phase, can be estimated by the minimum value at the end of smooth 
hybrid inflation, ~(v 2 / fJ 2 )<J c , where a c is the value of a at the end of smooth hybrid inflation. 
In our case, er c ~ 0.14. This gives the estimation of mean initial value ini at the onset of 
the new inflation 

v 3 

0ini ~ ~^3°" c - ( 16 ) 

With the total energy density at the reheating p re h, and assuming that ordinary thermal 
history after reheating, we can estimate the total e-foldings of inflation after the horizon 
crossing of the present Hubble scale: 

i 



67+ 2 ln JL + Ii n ^h. (17) 
3 Mq 3 Mq y ! 



New inflation 



During the new inflation, the dynamics of is controlled by 

^ ~ -C N v*<f> - 2^ngv n - 1 + l^ng 1 ^ 1 . (18) 



For v -C /i and /i given by (fTTj) . the number of e-foldings of the new inflation is 
estimated from Eq.t flTI) 1 as < 50. Assuming n = 4, we choose n = 4, Cn = 0.04, v = 
5.0 x 10 -4 , g = 2.0 x 10 -5 according to the dynamics of 0, in order to give Nn < 50. 

In this new inflation, the scalar potential has a negative minimum: 



V2 ( j^j " , V Nmin ~ -6 (-^ ) ^VLr (19) 



1 In principle, we must determine the decay rate of <fi into other light particles, including the standard 
model particles. Here, since it is enough to have rough estimation, we assume gravitationally suppressed 
interactions, as was done in 13 1. 
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This gives a negative cosmological constant after inflation. If we assume that there is another 
sector which breaks supersymmetry and that this negative energy density is canceled by a 
positive contribution from supersymmetry breaking, the scale v is related to the gravitino 
mass m 3 / 2 as 

m3/2 ~ ^tt (y)"'" 2 ' (20) 

Note that the present model v ~ 10" 4 gives unacceptably large gravitino mass, m 3 / 2 = 
1.2 x 10 8 TeV. However, we can abandon the relation between inflation scale v and the 
gravitino mass m 3 / 2 by assuming that the negative potential energy Vjvmin is cancelled not 
only by the contribution from SUSY breaking, but a constant term in the superpotential of 
another sector. In this paper, we assume that this is the case. 

We can estimate analytically the amplitude of primordial curvature perturbation for the 
mode crossing the horizon at the onset of the new inflation as 

K = — -= = — ■= a . (21 j 

2v / 37r0 ini 2V3n v 

This is larger than the amplitude at larger scales. Note that the slow-roll parameter 77^ is 
given by 

VN = ~C N - 2^n{n - l)g^- ^ ~C N . (22) 

Therefore, for the cosmo logically relevant scales which cross the horizon during the new 
inflation, the spectrum of curvature perturbation has an almost constant spectral index 
n s = 1 - 2C N < 1. 



III. NUMERICAL CALCULATIONS 

We now solve numerically the evolution of fluctuations of the scalar fields a, ip and 
until the end of new inflation, and calculate the spectrum of curvature perturbation after 
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the inflation. We use scalar potential approximated as 

V = V H + V N + Vj, (23) 
/ , ib 2m \ 2 ( a 4 ?b 2 ~ 

VH = l-H 2 + TZTTJ^-T, 1 + IT + 



4m^f2(m-l) J \ 8 2 

mVV / V V (m_1) m</> 2 "V 2 / 2 

+ 4 2m-l [jM J + 2 2m-lj^2(m-l) + ^m^m-l) 

V N = (^-^«) 2 -^V 2 , (25) 

/ 2 ^ \V / 2 \ , 

^ = + 4mM 2 (m -l) J Y~\ fi+ 4m M 2(m-l) ) V ^ ( 26 ) 

with the following parameter set: 

H = 2.04 x 10~ 3 , M = 1.17, v = 4.7 x 10~ 4 
m = 2 , n = 4, CV = 0.04. 



We adopt the linear perturbation formalism presented in 20j. We solve the evolution of 



perturbation in longitudinal gauge, in which perturbed metric is given by 

ds 2 = -(1 + 2<& A )dt 2 + a 2 (l + 2$ H )8 ij dx i dx j . (27) 

Evolution equations are given by 

k 2 d 2 V ( dV \ 

(Sep)'- + {3H + ri)(M; + -oScpi + V - - 5(fj - 2— + ^r,; $> H + = o, 

a 2 ^ dtpidipj V °<Pi J 

(28) 



(29) 



where roman subscripts run over 1,2,3, and <pi, (p 2 , v?3 stand for a, ip, 4>, respectively. 5fi 
means fluctuation of corresponding scalar fields. 

We calculate curvature perturbation on uniform- density hypersurface (: 

C = + H-^h) +(l+ Q 2 f m P . ) 9 B , (30) 

3p + p \ 9a 2 H z p + pJ 

which is constant and equal to comoving curvature perturbation 1Z on superhorizon scales. 
Here, p and p are total energy density and pressure of background, respectively. We introduce 



decay of scalar fields into radiation, with homogeneous and time-independent decay rates in 



the manner introduced by 20j]. This is done by adding an extra friction terms to the equation 
of motion of scalar fields. We treat radiation as a thermal bath, which has contribution to 
energy density and pressure of background. The energy density of radiation p r decays in 
proportion to a~ 4 , while energy is injected from decay of scalar fields: 

^^^ = r^ + T^ + T^. (31) 

We set decay rates IV and to being negligibly small during the whole calculation : r CT = 
r^, = 10~ 10 M G <C H N <C H H . On the other hand, since gives reheating temperature, 
should be far smaller than this. In our calculations, we regard to be completely negligible: 

r = o. 

In Fig{I|-(a), we show the spectrum of primordial curvature perturbation calculated nu- 
merically. At ko = 0.002[Mpc~ 1 ], this spectrum has 

ft = 4.9 x 10" 5 , n s = 1.053, = -0.032. (32) 

a In A; 

This is not the best-fit value of the WMAP three year result, but within la range. 

On large scales, the spectrum has a desired shape and amplitude, although it does not 
agree with the best-fit parameters. The amplitude becomes large at k = l[Mpc ], and is 
connected to the spectrum produced by the new inflation, which has larger amplitude at 
k ~ lOfMpc -1 ]. According to large scale structure observations, the spectrum is constrained 
to be almost flat on scales larger than k ~ l[Mpc -1 ]. This is satisfied well by the spectrum 
in FigEEKa). 

A sequence of sharp peaks is seen between ko = 200[Mpc _1 ] and fco — 1500 [Mpc -1 ]. This 
originates from parametric resonance, as discussed in the next section. The height of the 
peaks is determined by competition between parametric resonance and decay of the scalar 
fields a, ip. When the decay rates are large, the fluctuations (= scalar particles) decay before 
the amplitudes are amplified through the parametric resonance, which results in suppression 
of the peaks in the power spectrum of the curvature perturbation. This is seen in FigJT]-(b) 
where the decay rates of a and ip are assumed to be 6 x 10~ 8 Mg. 

In other words, we design the height of the peaks by choosing appropriate couplings 
between a (ip) and the standard model particles. 

If we take observation of Ly-a into account, the amplitude of the perturbations should be 
small sufficiently for k < lOfMpc -1 ]. The spectrum shown in Figj2] meets this requirement. 
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10" 4 10" 3 10" 2 10" 1 1 10 10 2 10 3 10 4 



k [Mpc" 1 ] 

FIG. 1: (a) The spectrum of primordial curvature perturbation produced by smooth hybrid new 
inflation, calculated under the parameter choice (|32l) and negligibly small decay rate T = 1CP 10 Mg 
of inflaton fields a and if). The largest peak is located at k = 810[Mpc -1 ]. (b) The spectrum 
calculated under the same parameter choice for (a), but larger decay rate T = 6 x 10~ 8 Mg. 

It is given under the parameter set: 

H = 2.04 x 1(T 3 , M = 1.12, v = 5.0 x 10~ 4 

m = 2 , n = 4, C N = 0.04. (33) 

At ko = 0.002[Mpc _1 ], this spectrum has 

n = 4.9 x 10~ 5 , n s = 1.10, — ^- = -0.026, (34) 

am k 

with n s and dn s /d\nk at the edge of la range of the WMAP three year result. 

IV. PARAMETRIC RESONANCE IN SMOOTH HYBRID NEW INFLATION 
MODEL 

In this section, we will investigate the resonant amplification of fluctuations of the scalar 
fields 5a k and Stpk, an d consequent amplification of 54>k- This causes a strong peak in the 
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k [Mpc" 1 ] 

FIG. 2: The spectrum of primordial curvature perturbation which is sufficiently small for k < 
10[Mpc -1 ]. This spectrum is calculated under the parameter choice (|33p and negligibly small 
decay rate T = 10 _10 Mg of inflaton fields a and tp. 

spectrum of primordial curvature perturbation. Since these two regimes have very different 
time scales, we will describe them separately. 

A. Parametric resonance regime 

After the smooth hybrid inflation ends, a and ip begin to oscillate about their respective 
minima. We show these oscillating backgrounds in FigJSJ Hereafter, we normalize the scale 
factor to be a = 1 at present. 

The evolution is very complicated. Because the cross terms which includes both a and 
if) exist in the scalar potential, Sa^ and 5ipk contributes to the evolution of each other. 
Moreover, at the beginning of the oscillatory phase, the second order terms in oscillating 
backgrounds o and if) contribute to the evolution as well as the first order terms, which 
makes analytical understanding of parametric resonance difficult. 

Retaining only the first order terms, and neglecting metric perturbations, evolution equa- 
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200 400 600 800 1000 
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FIG. 3: Oscillation of backgrounds a and ip. Time variable t is indicated in unit of m~ , here 
m = m CT ~ m^p. t = is taken to be In a = —114.5, at the time shortly before the oscillatory phase 
begins. 



tions of 5<Jk and dipk are given by 
(5a k y + 3H(Sa k y + 

(W + 3#(W + 



"fc 2 


+ m 2 a + 


24/i 2 


a 2 


M 2 


p 2 
a 2 


+ ml + 


36/i 2 
M 2 



5(Jfc + 



24/i 2 
M 2 



^/JlMadipk = 0, 



(35) 



6fa + ^-y/jMa6a k = 0, (36) 



and retaining lowest order, evolution equations of backgrounds cr and ^ are given by 



a + 3H& + mia = 0, 



(37) 
(38) 



ijj + 3#^ + ml4> = 0. 



In order to understand their behaviors approximately, we will recast Eqs. fl35l) and ([36 
into following form: 



5a" + [A a - 2q U(T cos (2(1 + d)z + A)] 5a - 2q^ cos{2z)5ip = 0, 
5a" + [Af - 2q H , cos (2(1 + d)z + A)] 5ijj - 2q a1 p cos(2z)5a = 0. 



(39) 
(40) 
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Here, the prime represents the derivative with respect to the variable z defined as 2z = 
m a t — 7r/2. A is a possible phase, and 1 + d = m^/m^, which satisfies d <C 1. Furthermore, 
we approximated the oscillating background as 



a ~ — T,sm(m a t), (41) 
ip ~ —ty sm(m^t + [phase difference]). (42) 

We also rescaled the perturbations ba k and 5?pk as 5a = c?l 2 bo k and 5ip = a 3 ^ 2 5ip k , and 
neglected 0(H 2 ) terms. The coefficients A a , A^, q ua , q a1 p, q^ a , and q^ are given by 

4k 2 ( k 2 \ 

A a = A + — — , ^ = 4 + 4 — — + d , (43) 
a z mi \ a z mi ' 



4V^M 4^M T^/JIM 

Qct<7 = A/f2 2 W, q&tp = „ 2 £, = „ 2 W. (44) 
M 2 m% M 2 mi M 2 mi 

a a a 

At the beginning of the oscillatory phase, amplitudes are estimated as E ~ $ ~ 0.01. 
Neglecting d < 1, we get A a ~ ~ 4-^r + 4 > 4, g CT(T ~ q ai) ~ 0.62, and ~ 0.93. 



Eqs.(|HS]) and (SO} look like the Mathieu equation |2J] 



x" + [A-2q cos(2z)}x = 0. (45) 

Thus, we expect instability solutions. 

To confirm the existence of the instability, we will show some results of numerical calcu- 
lations. First, we have solved the coupled evolution equations (I35p - (I38p numerically. In 
order to trace evolutions of perturbations 5a k and 5ip k and to see the occurrence of resonant 
amplification, we neglected expansion of the universe: we put H = 0. 

Figure H] and Figure [5] show time evolutions of power spectra 2 of scalar field perturbation 
V S a k and V s ^ k under Eqs. (EBJ - (JMD- We use S = 1.08 x 10~ 2 and = 1.16 x 10~ 2 , which are 
evaluated at In a = —114 from the full numerical calculation of smooth hybrid new inflation 
given in Sec. III. Initial amplitudes \5a k \ and \5ip k \ are set to give \ogVs ak = and logVgf k = 
0, respectively. We concentrate on two modes: k/a ~ 3.5 x lO -1 ?^ and k/a ~ 7.0 x 10~ 7 m a . 
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For Fourier modes <?£ of a perturbation g{x), the power spectrum V g is defined by 
Here, (• • • ) represents an ensemble average. 
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These modes correspond to k — 1000 [Mpc -1 ] and k = 0.002 [Mpc -1 ] at present. We can see 
that resonant amplification takes place. Rate of exponential amplification varies depending 
on k. We also find that 5a^ and 5ipk show almost identical evolution. 

Figure [6] shows scale dependence of amplification under the same configuration as that 
used for FigJH We indicate only the amplification of Sa^, at the time t = SOOm^ 1 . Since 
5ak is oscillating rapidly, we use Vs<r k | pea k, the peak value of oscillating Vs ak around at some 
specific moment, which is approximately equal to the amplitude of oscillation. We can 
see that small momentum modes k 2 /a 2 <C m 2 show identical amplification. Most efficient 
amplification takes place at k/a ~ 3.5 x lO -1 ?^. 

Next, in order to compare the above discussion with the actual time evolution of smooth 
hybrid new inflation model, we numerically calculate the evolution of the fluctuation 5<7k, 
in the same way as described in Sec. Ill, solving all linear evolution equations of scalar fields 
Eq. fTSBl) and metric perturbation Eq. (l29]) under the scalar potential V given by Eqs. fTSBl - 
(f26|) . in the expanding universe. 

Figure [7] shows the mode dependence of efficiency at various times while the resonant 
amplification takes place. We use Vs<r k | P cak, for the same reason that we did so for FigJBl 
Since all modes of fluctuations of scalar fields decrease because of the Hubble friction term, 
the amplification is suppressed during the expansion of the universe, so that modes around 
the strong peak are amplified significantly in contrast to other modes. At £ = 15m" 1 , we 
can see efficient amplification of low-momentum modes. Since the effective potential of a is 
tachyonic at the end of hybrid inflation, all modes with k 2 /a 2 <C m 2 , are amplified efficiently 



by tachyonic instability 22, 
perturbation at this time. 



231 ]. No significant peak appears in the spectrum of scalar 



B. Forced oscillation of 5(pk 

Now, we turn to the evolution of 5<pk- We will focus on subhorizon modes at the beginning 
of oscillatory phase: k/a > = /i 2 /v^3 with In a = —114.5. For superhorizon modes, 
curvature perturbation is already frozen out, and evolution of 84>k after smooth hybrid 
inflation is irrelevant. 

During the smooth hybrid inflation, evolutions of Sa^jSip^ and 5<pk are all dominated 
by k 2 /a 2 term. They decay like |<5cr| 2 oc a~ 2 from the same initial condition given out of 
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t [ma -1 ] 

FIG. 4: Time evolutions of Vsa k under Eqs. ([35|) - ([38]) . for k/a ~ 3.5 x 10~ 1 m a and k/a ~ 
7.0 x 10- 7 m a . 
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FIG. 5: Time evolutions of Vgip h under Eqs.(|35p - (|38p . for k/a ~ 3.5 x 10 1 m CT and k/a ~ 
7.0 x lQ- 7 m a . 
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FIG. 6: Scale dependence of amplification under the same configuration as that used for Figf5]and 
FigEJ We show /c-dependence of (Vs(r k Ipeak) / (Vs<r k \peak,t=o)- Most efficient amplification takes 
place at k/a ~ 3.5 x lQ~ l m a , at the time t = 500m" 1 . 



quantum fluctuation. So, we can estimate that Sa k , 8ip k and 
at the beginning of the oscillatory phase. 

Evolution equations for Fourier modes 8cj) k are given by 



have the same amplitude 



(<ty fc )" + 3H(5<p k y + 



k 2 4/i 3 , 2 

V 2 + ^ 



b h - 2 ~g^Mv 2 (iP8a k + a5^ k ) = 0. (46) 



No resonant amplification due to oscillating a and ip occurs, since (4/z 3 /M)\I/ 2 <C k 2 /a 2 <C 



m 



Let us assume sufficiently large amplification of 5a k or 5ip k occurs so that (ip5a k + <j5ipk 
terms dominates Eq.( l4*6|) . The condition for that is given by 

\k 2 



< 



^^Mv 2 m5a k 



(47) 



In our case, this requires more than 10 4 amplification of 5a k relative to 5<f) k . Since we have 
more than 10 5 amplification at the resonant peak, this condition is satisfied very well. 

Once the (ip5a k + aSipk) term dominates Eq. (|46|) . 54> k undergoes forced oscillation with 
the source term (ip5a k + a5ifj k ). One can expect that this results in an amplification 
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FIG. 7: Scale-dependence of the resonant amplification. We show /c-dependence of 
a k o 3 1 peak, t=o ) • The bottom line indicates t — Ihm^ 1 , the rest show various t 
from t = 100m" 1 (next-to-bottom line) to t = lOOOm" 1 (top line), varying by lOOm^ 1 between 
each lines. At t = 15m" 1 , amplification by tachyonic instability can be seen: all modes k/a -C m 
are equally amplified. After that, resonant peak appears around at k ~ 1000[Mpc -1 ]. 



of S(pk, but there is another subtlety. In our case, ip and 5ipk oscillate with frequency 
about m^,, while a and 5au oscillate with about m a . Since these two frequencies are 
slightly different from each other, tpSak and aSipk show log-period oscillatory behavior as 
(C/2) [cos(m a St) — cos(m CT (2 + 5)t)}, where C is the amplitude. Then, S(pk has a solution 
like 

Ca 2 

5<fik ~ —^■cos(m a 5t). (48) 

As a result, S(pk shows long-period oscillation, whose frequency is estimated by Am. This 
induces large oscillation of the amplitude of 5<pki while the mean value of oscillation is still 
decaying. Meanwhile, the source term decays with I^Okl oc a -3 , since a and ip behave as 
massive scalar fields. On the other hand, H takes a constant value after the beginning of the 
new inflation. Therefore, contribution from Hubble friction term eventually becomes relevant 
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and 8<pk ceases oscillation. Afterwards, 5<pk decays until the horizon crossing, and then it is 
fixed. 5<fik at that time determines the amplitude of primordial curvature perturbation: 

Vn= {j) ^ (49) 
Since the amplitude of 5<frk is determined by the phase of large oscillation, /c-dependence of 
54>k at horizon crossing is oscillatory. Thus, resultant spectrum of the primordial curvature 
perturbation has a strong peak, "sliced out" from the single resonant peak generated by 
parametric resonance of a and if) as shown in FigJHl where we show the time evolution of 
5(f>k in this regime. 

-4r 
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-6 \ 

-7 | 

-e- l 
to 

^ -8 ! 

O r 

bO o | 

r2 S 

-10 ! 

-li ! 

-12 L 1 

-13 ^ ! 

2000 4000 6000 8000 10000 12000 

t [mo" 1 ] 

FIG. 8: Long-period oscillation of Vs<j> k for the mode of resonant peak: k = 810[Mpc -1 ]. Sharp 
valley indicates that the oscillation crosses zero. The period is estimated to be At = 2tt/ Am ~ 
2.6 x 10 3 m~ 1 , which agrees well. 

Let us reexamine the mechanism which induces this resonant peak by comparing estima- 
tion of the strength of peak and the result of numerical calculation. 

At the beginning of long-period oscillation, the equation of motion (1461) is dominated 
by (k 2 /a 2 )S(pk term and (ipSa^ + crSipk) term. Therefore we can estimate the amplitude 
of long-period oscillation of 8cf>k from above discussion. The first peak of V^ k is at about 
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t ps 1800m" 1 , lnasi —111.8. At that time, the amplitude of 8<Jk is given by Vg ak ~ 4.5 x 10 2 , 
while the amplitude of oscillating background is \I/ ~ 3.2 x 10 -4 , according to the numerical 
calculation. We can estimate the amplitude of the long-period oscillation of 5(j)k, in terms 
of the power spectrum for k = 810[Mpc~ ], 

a 4 4uV# 2 

v*. - hhr~ Vsai s 17 x 10 (50) 

Hereafter, we will concentrate on amplitudes 8(Tf. and Scfik of the oscillations of 5cr& and 50^, 
respectively. We approximated that the a5ipk term gives the same contribution as that of 
i/j5<Jk term. The numerical calculation gives V s j >k ~ 9.3 x 10 5 . Our estimation is very rough, 
but gives the same order of the result of numerical calculation. 

At the time the long-period oscillation ceases, this estimation underestimates the am- 
plitude V 5 ^ k . At t ~ 750077^ 1 and In a ps —108.8, the long-period oscillation ceases, with 
V S a k ~ 3.7 x 10" 2 and * ps 3.5 x 10~ 6 . These give V s i ~ 9.9 x 10" 8 , while the numeri- 
cal result is V s i k ~ 1.0 x 10~ 6 . This is because of the complexity of the evolution: since 
k/aH ~ 0(1), the evolution is in the transition regime between sub horizon to superhorizon, 
and hence the evolution of \8<f>k\ marginally freezes. 

On the other hand, we can verify that 8(pk at the horizon crossing determines the curvature 
perturbation. At the horizon crossing, numerical calculation gives V S( j >k ~ 3.0 x 10~ 7 . At 
that time, H ~ 1.3 x 10~ 7 and <fi ~ —1.3 x 10~ n . We can get 

2 

I 

Ipeak~ I 1 I <5</>fclcross 



^1 ) V si \^29. (51) 



This gives log TV < 1.5 at the peak, which is the same order of the result of numerical 
calculation log 7^ ~ 1.15. 



V. PBH FORMATION FROM RESONANT PEAK 

The presence of the large sharp peak around k ~ 10 3 [Mpc _1 ] may lead to interesting 
consequences, one of which is the formation of primordial black holes. From the peak 
position k = 810[Mpc _1 ] of the spectrum shown in FigJT]-(a), the mean mass of resultant 
PBHs is estimated as 

m BH = 1 x 1O 9 M . (52) 
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Here we assume that the whole mass within the overdense region when it enters the horizon 
collapses into a black hole. The abundance of such a massive PBH is constrained by assuming 
that they must not overdose the universe {24)]. In this case, initial mass fraction (3 is 
constrained to be (3 < 4.6 x 10~ 6 . 

We estimate the black hole abundance, following [25] based on numerical simulation of 



PBH formation from an overdense region 261 ] . 



Assuming a spherically symmetric overdense region, two conditions must be satisfied 
for the formation of PBHs from the overdense region. First, the overdense region must 
have larger density than a critical value. Under the linear approximation, this is given by 
curvature perturbation at the horizon crossing 7Z(t, r) as 

K(t,0) > 0.67 : CI. (53) 

The other independent condition is given by 

- 0.72 < X(r) < -0.28 : C2, (54) 

where X{r) = at the horizon crossing. The latter condition requires that the excess 
mass around the overdense regin is sufficiently large. The initial mass fraction f3 of PBHs 
can be identified with the probability P(C1 PI C2). Assuming that the density perturbation 
has a Gaussian probability distribution, the probability P(C1) is given by 

P(C1) = f — L_ exp dS, (55) 

where a-ji is the variance of TZ at the horizon crossing. For the spectrum shown in FigiU-(a), 
we can get an = 1.9, which results in (3 ~ 0.35. 

On the other hand, under the presence of a peak with lZ(t, 0) = 0.67, the conditional 
probability of (C2) is found to be 0.5. Consequently, the initial mass fraction of PBHs is 
estimated to be 

/3~P(ClnC2) -0.18. (56) 

This is far larger than the constraint f3 < 4.6 x 10~ 6 . 

If we introduce larger decay rates of inflaton fields a and ip, resultant PBH abundance 
can be acceptably small. The spectrum shown in FigJT]-(b) gives (3 ~ 5.6 x 10~ 9 , which is 
below the bound (3 < 4.6 x 10~ 6 . 
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VI. CONCLUSION AND DISCUSSION 



We have reexamined the smooth hybrid new inflation model, which was designed to 
reproduce the running spectral index suggested by the WMAP three year result, by numer- 
ical calculation of the perturbation. We have confirmed that this model can reproduce the 
spectrum within la range of the WMAP three year result. 

In addition, we find strong peaks on smaller scales, which originate from the amplification 
of perturbation 5(f) and 5if) by parametric resonance. Since there are interactions between 
a, if) and <f>, this amplified perturbation is transfered to 5(f). 5<f> begins large-amplitude and 
long-period oscillation, and freezes out at the horizon crossing, determining the primordial 
curvature perturbation. The resultant curvature perturbation depends on the phase of 
oscillating 5<f>, so that the resonant peak consists of several steep peaks and valleys. The 
strength of the resultant peak can be controlled by the decay rate of the first inflaton, and 
the comoving scale of the peak can be controlled by the scale factor at the beginning of 
oscillatory phase. 

Such a peak in the spectrum of primordial curvature perturbation can be a source of 
PBHs. Due to the strong peak, mass of generated PBHs can be estimated by the horizon 
mass at the horizon crossing of the scale which corresponds to the resonant peak. 
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